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We study the p meson in a uniform magnetic field eB using a holographic QCD-model, more 
specifically a D4/D8/D8 brane setup in the confinement phase at zero temperature with two 
quenched flavours. The parameters of the model are fixed by matching to corresponding dual field 
theory parameters at zero magnetic field. We show that the up- and down-flavour branes respond dif- 
ferently to the presence of the magnetic field in the dual QCD-like theory, as expected because of the 
different electromagnetic charge carried by up- and down-quark. We discuss how to recover the Lan- 
dau levels, indicating an instability of the QCD vacuum at eB = m^ towards a phase where charged 
p mesons are condensed, as predicted by Chemodub using effective QCD-models. We improve on 
these existing effective QCD-model analyses by also taking into account the chiral magnetic catalysis 
effect, which tells us that the constituent quark masses rise with eB. This turns out to increase the 
value of the critical magnetic field for the onset of p meson condensation to eB k, 1.1 m^ k, 0.67 GeV^. 
We briefly discuss the influence of pions, which turn out to be irrelevant for the condensation in the 
approximation made. 



I. INTRODUCTION 



Recently the interest in magnetically induced QCD effects has grown, as a magnetic field offers a controllable 
parameter that gives rise to new unexpected results that might lead to novel insights in QCD, e.g. the chiral 
magnetic effect, which is related to CP-violating processes chiral magnetic spiral /wave, a split between the 
deconfinement and chiral symmetry restoration temperature, etc. Since at the LHC very strong magnetic fields 
(of the order of 10^^ Tesla) are expected to arise in heavy ion collisions |l2]-0, this creates the perfect setting for 
experimental searches for possible interesting QCD effects in a strong magnetic field. This has resulted in a lot of 
activity in this relatively young research field. We refer to [7-59 1 for a selection of relevant literature. An interesting 
review is Il60l . During the cosmological electro weak phase transition, gigantic magnetic fields should have occurred 
as well Il63ll . 

Although realistic magnetic fields might have a short lifetime and/or be spacetime dependent, we will consider 
them fixed here, as is usually done. 

One particular new QCD phenomenon, put forward by Chemodub in ||30l |3T| and further explored in 1 32^341 |36ll 
is the possible instability^ of the QCD vacuum towards condensation of charged p mesons when a sufficiently 
strong magnetic field is present, eBc ~ m^, somewhat analogous to a possible W^-boson condensation in the 
electroweak model Il35ll64ll65ll . In [30 31 [ a phenomenological approach was adopted to argue the existence of a 
p meson condensate, respectively using an effective quantum electrodynamics Lagrangian for the p mesons and a 
NJL-model. Lattice evidence in favour of this effect appeared in [40 J. A consequence of such a condensation would 
be that the QCD vacuum at zero temperature becomes superconducting at sufficiently large magnetic field, forming 
a quite exotic phase. In addition, it should be a kind of anisotropic superconductor, see [30-32| for more details 
on this. Quite recently, a similar instability was studied in a holographic toy model [55 56J. Other examples of 
magnetic instabilities of bulk charged vectors in an AdS/CFT context can be found in [|61H62ll . 

We set out to give further evidence for a magnetically induced p meson condensation using a holographic approach. 
Some preliminary and incomplete results were already presented in [50 1. The investigation of the condensation 
itself and ensuing conductivity properties of the vacuum will be relegated to future work. Here we intend to take 
the first step in this program, and show that a condensation should occur as we encounter a tachyonic instability in 
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the p meson sector. We shall rely on the much studied Sakai-Sugimoto model Il66ll67ll . 

Holographic QCD-models give a description of hadronic physics through a dual supergravity theory in a higher- 
dimensional world. This working philosophy is based on the original Maldacena conjecture |69|. The duality 
is valid for a large number of colours and at large 't Hooft coupling A where QCD itself is unmanageable, thus 
providing an analytical setting for studying non-perturbative QCD effects. The Sakai-Sugimoto model in particular 
Il66ll67l uses a D4/D8/D8 brane setup which manages to reproduce much of the low-energy physics of quenched 
QCD in the chiral limit, such as confinement and dynamical chiral symmetry breaking at low temperature, vector 



meson dominance, pion effective action, etc. A short review of the model is presented in Section II A We will 
discuss the effect of the magnetic field on the p meson mass in both the original antipodal embedding considered 
by Sakai and Sugimoto (where the flavour branes merge at a value uq of the holographic radius equal to the cut-off 
of space uk) and the more general non-antipodal embedding (mq > uk)- The latter enables us to take into account 
constituents of the p meson, at least to some level. We extend the work of Sakai and Sugimoto on the numerical 



fixing of the holographic parameters for the uq — Uj^ case to the case where uq > Uj^ in Section II B In order to get 
results in physical units we have fixed the free remaining parameters in the holographic model by matching them to 
phenomenological or experimental values for the constituent quark mass, the pion decay constant and the p meson 
mass nip, in absence of a magnetic field. Let us mention here already that throughout this article the notation mp 
will be solely used to denote the p meson mass at zero magnetic field. 

In the Sakai-Sugimoto model at zero background magnetic field, the flavour D8- and D8-branes form a stack of coin- 



ciding branes in the D4-brane background. In Section II C we summarize how to introduce a background magnetic 



field eB into the model. Only in the mq > case there is an effect of eB on the embedding. Considering two flavours, 
we allow for the possibility that the up- and down-flavour branes respond to the magnetic field in a different way, 
which they turn out to do. This simply reflects in a geometrical way the up- and down-quarks coupling to the mag- 
netic field with different electric charges. In other works like ll47ll49ll5TH53l , the influence of a magnetic field on some 
features of the Sakai-Sugimoto model was already investigated, however mostly for the = 1 case. The chiral 
magnetic effect was studied for example in |49,.,53J, whereas a part of the phase diagram in terms of a magnetic field 
and isospin chemical potential in [51 1. Since we are interested in charged p mesons, we need multiple flavours in our 



study, and thus have to extend on existing Nf = 1 literature^. We find in Section II D among other results, a heavier 
constituent mass for the up-quarks than for the down-quarks in the presence of a magnetic field (in the uq > Uk case). 



In Section III A we derive the mass equation for the p meson, recovering the Landau levels of fSD] that indicate the 
instability of the vacuum, as the effective masses of the charged p mesons' polarizations that are transverse w.r.t. the 
applied magnetic field become imaginary in the lowest Landau level when the magnetic field reaches the critical 
value eBc = m^. An obvious drawback of the phenomenological model used in fSO] is the fact that the p mesons are 
regarded as point particles without internal structure. The magnetic field is so strong that it is not evident to trust 
that the underlying quark structure will be of no importance to the story. Using a NJL model allows to consider 
internal quark dynamics IISTI but the ensuing p meson dynamics is then perhaps somewhat less phenomenologically 
valid than that of [30] . 



In Section III B we take one effect of the constituents of the p meson into account through the effect of chiral magnetic 
catalysis, which tells us that the constituent quarks get heavier in the presence of a background magnetic field, as 
will the p meson itself (apart from the Landau shift lowering the effective mass). The resulting effective p meson 
action resembles that of [30] but with an effective eB-dependent mass due to internal quark effects. This leads 
to a slightly larger value for the critical magnetic field, eBc ~ l.lmi, at the onset of the predicted condensation. 



This increase is in qualitative agreement with the lattice data of [40J. In Section IIIC we discuss that pions, in a 
leading approximation in the inverse 't Hooft coupling, should not disturb the charged p meson condensation, as 
the magnetic field induces a coupling only between the pions and non-transverse, w.rt. the magnetic field, p meson 
polarizations. We summarize in Section llV] and point out how to improve our analysis in future work. In the 
Appendix, we have collected some technical details. 



^ In 1521 the two-flavour Sakai-Sugimoto model was used. The embedding however was the same for up- and down-flavour branes, even in the 
presence of a magnetic background, because only the Uq = case was considered there. 
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II. HOLOGRAPHIC SETUP 



A. Review of the Sakai-Sugimoto model 



The Sakai-Sugimoto model 
brane background 



1671 involves a system of Nr pairs of D8 — D8 flavour probe branes placed in a D4- 
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where dO^, and V4 = 87t^/3 are, respectively, the line element, the volume form and the volume of a unit four- 
sphere, while R is a constant parameter related to the string coupling constant gs, the number of colours Nc and the 
string length £s through R^ = nggNci^- This background has a natural cut-off at u = and the QCD-like theory is 
said to "live" at m ^ 00. Imposing a smooth cut-off of space at m = Uk determines in a unique way the period of 
t: 
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with Mk the inverse radius of the T-circle. In this paper we will work with 3 colours, Nc = 3, and 2 flavours, = 2. 
We stress here that we ignore the back reaction of the flavour branes on the background, an assumption which is 
in principle only valid for Nc ^ Njr, which means that we are working in a holographic analogue of the quenched 
approximation. Unquenching the Sakai-Sugimoto model is a difficult task, see IItTI . 



U/N 





FIG. 1: The Sakai-Sugimoto model: antipodal {hq = ug^) and non-antipodal (mq > mx) embedding. 



On the stack of Njr coinciding D8 — D8 flavour pairs, there lives a U{Njr)i x Lf (N^)r gauge theory, which is inter- 
preted as corresponding to the global chiral symmetry in the dual QCD-like theory. The cigar-shape of the (m, t) 
subspace of the D4-brane background enforces a U-shaped embedding of the flavour branes, encoded in m(t). This 
embedding represents the dynamical breaking of chiral symmetry U{Njr)i x U{Njr)R U{Nj-) as the merging of 

the D8-branes and D8-branes at m = uq. The asymptotic separation L (at m 00) between D8- and D8-branes, 
indicated in Figure [l| is related to uq as 

^ [u' = du/dr) 
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In the original setup of Il66ll67ll the flavour branes merge at the tip of the cigar, mq = "K/ the embedding is antipodal. 
The more general non-antipodal embedding with uq > Uf^ models non-zero constituent quark mass |70| 

IllK Vf{u) 
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with oil = the string tension. In this setup, unlike in the mq = case, it is possible to model chiral magnetic 
catalysis Il72l : a magnetic field boosts the chiral symmetry breaking and hence the constituent quark masses. 
This phenomenon is by now generally accepted to occur in QCD, having received numerous confirmations IITTl - 
IT5llT9ll37l l4?-^,'47l l48ll73ll , albeit that the most realistic QCD lattice data display a richer behaviour than originally 
anticipated [44, 45], probably related to quark mass/pion effects^. For the Nj = 1 case of the Sakai-Sugimoto 
model, the catalysis was already observed in [ [47| . The current quark masses are always zero in the Sakai-Sugimoto 
model, meaning that we are working in the chiral limit. This is closely related to the absence of a direct holographic 
description of the chiral condensate in the Sakai-Sugimoto model, which can only be overcome at the cost of great 
complication Il74l . 

The U{Nj:) gauge field Am (x'', «) {m = 0, 1, 2, 3, m) living on the D8-branes describes mesons in the boundary field 
theory. The action for this gauge field is given by the non-Abelian DBI-action^ ll75H78l 



Sdbi = -Ts J d^x 2 du J 64 e-1> STr^- det [g^?^ + {2n a' (5) 

where Tg = l/((27r)^£g) is the D8-brane tension, the factor 2 in front of the m -integration makes sure that we integrate 
over both halves of the U-shaped D8-branes, STr is the symmetrized trace which is defined as 

1 

STr(Fi ■ ■ ■ F„) = — |^Tr(Fi ■ ■ ■ F„ -|- all permutations), (6) 

g^ifi the induced metric on the D8-branes, and Fmn = dmA„ — d„Am + [Am,A„] = F^^^jt" the field strength with 
anti-Hermitian generators 

f = -'-il,crt,a2,a3\ Tr(i''i^) = [t\t''] = eabc^. (7) 
The parameters R,gs, £s, Mk, mk and 't Hooft coupling A = ^y^Nc are related through the following equations: 

Without loss of generality^ one can moreover set |M|^£s = \ II67I which is the same as saying 



Consequently 



Following ||67|| , we define 



for notational convenience. 
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^ In particular, regions of magnetic "inhibition" at finite temperature were reported, with analytical insights discussed in f24M25!. 

* The full non-Abelian generalization of the DBl-action is not known. The STr-prescription by Tseytlin which we will use is valid up to fourth 
order in the field strength, with deviations starting to appear at order 17911801 . The DBl action ignores derivative terms including [F, F] ~ 
[D, D]F terms (ambiguity of the slow-varying fields approximation in non-Abelian case) |75|. Because the 't Hooft coupling A is large in the 
validity regime of the sugra/QFT duality, we ignore in a first approximation the Chern-Simons part of the action in the analysis, being a factor 
1/A smaller than the DBI-action |67, 90 [. We will come back to the discussion of this part of the action later. 

^ All physical results are independent of the choice of Af^ . 
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B. Numerical fixing of the holographic parameters 

In the model with antipodal embedding there are six parameters, R, k, £s, Mk, gs and L, related to each other through 
the four relations l|8|-d9l. In [67] the remaining independent parameters and k were fixed to GeV units by 
matching to the QCD input values 

/„ = 0.093 GeV and nip = 0.776 GeV, (12) 

resulting in 

Mk w 0.949 GeV and k = w 0.00745. (13) 

21671'^ 

We extend this to the model with non-antipodal embedding where there are seven parameters, R, k, is, Mk, uq, gs 
and L. These are completely determined by the four relations l|8|-l|9}, and the three extra requirements that the 
computable numerical values for the constituent quark mass niq, the pion decay constant and the p meson mass 
ntp, in absence of magnetic field, match to the phenomenologically or experimentally obtained QCD input values 

= 0.310 GeV, = 0.093 GeV and mp = 0.776 GeV. (14) 

Adapting the analysis of Sakai and Sugunoto ||66l to the more general case uq > uk, see also 1821 |83l, we derive the 
mass eigenvalue equation for the vector meson sector and the expression for f„ as functions of Mk, mq arid ^, and 
thus indirectly (through the relation ||3}) as functions of the three unknown parameters Mk, L and k. 

In the confinement phase, the non-Abelian DBI-action ||5j becomes, to second order in the field strength, 

Sdbi = V jd^xl du Tr {u-'/^j'/^R%,F>''' + lu^'S~^'^n''"hu^vu} + 0{F% (15) 



with 



and 



V = TsVig-'R'/^l{27z^')^ 



(16) 



and where the STr is replaced by Tr as 



STrift'') = Trift''). (18) 
Assuming the flavour gauge field components A,, (x'', m) can be expanded in complete sets {ij)„ («)}„>!, 

A^{x,u) = J2 Bl;'\x)ip„{ii) (19) 



n>l 



withF^™' = d^tSi'"'' - dyBl'"\ FftuFvu = B^!"^ bI"^ dufmdu^n 



the part of the action ijisjl in b|,"^ (x) reduces to the effective 4-dimensional action 

S = - I |1f/,;:'''f(")''™ + ^m2Bi")"B(«)'"'| , (20) 

describing vector mesons bJ"' with masses ifiu/ if the ipni^u^ are subject to; 

V / du U-^'W^R^,n^n = -5„,n, (21) 
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and 

2V duu^'^j-^'^{du^,n)[Mn) = ^mlSmn {m,n>l). (22) 

The conditions | [2T| and | [22| combine, using partial integration (and hereby implicitly assuming the i/',, to be normal- 
izable fimctions)7to the eigenvalue equation 

- u^'^^-^'^d,{u^'^j-^'^du^„) = R^inl^. (23) 

To include pions, following the original discussion in |66|, it is convenient to introduce a new coordinate z that is 
related to u through 

= ul + uqz^, (24) 
going from — oo to +oo along the flavour branes and thus allowing the description of both boundaries at m ^ oo of 



the U-shaped flavour branes. In this new coordinate z, the action 1 15 1 reads (denoting u{z) as Uz for readability) 

//■oo (o ^ly^/-^ 1 

/ dz Tr <^ Uoul/WR%.vFf'' + - ^V^'FhzFuz + 0{F% (25) 
J— oo \ J uq 'y I 

with 



and the condition d22 



^ |z I / z 

^ y Mf(wi - m|) - (W^ - W^w|)' ^^^^ 



V 3 u^^^ 1 

- / dz '—{dzip„,){dzip„) = ^m\8,nn {m,n > 1). (27) 

2 j-co Mo 7 2 

The flavour gauge field component^ Az is expanded in the complete set {(pn{z)}„yQ'- 

A^(x,z) = ^bJ"'(x)!/^„(z), Az{x,z)=J2'P^"^i^)'Pniz) (28) 

tt>l n>0 

^ F^,F/"' = FpF(")f'iPmip„ , F^zFyz = (9,,</'('"V<" " Blr^dz^,„){d,cp^"^cl,n - Bi"^dziPn)- 

Demanding a canonical normalization of the kinetic term for the <^(") (x^) fields in the effective 4-dimensional action 
then leads to the orthonormality condition 



dz —ui'^i-'(Prn(pn = iz&mn {m,n > 0). (29) 
^ J-oo Uq Z 



From the last condition | |29| and | |27| it follows that 

^„ = m^^dzipn for M > 1 (30) 



and since -L (^z^n for all n > 1, 



we can set 



dz — M^/^y->o3z!/'n = 0, (31) 

Z J-oo Mo 



' The gauge field components along the four-sphere are assumed to vanish, = 0, and Af, and Az are assumed to be independent of the 
four-sphere coordinates. 
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with the normalization constant c determined by 



^c^ r dz = ^. (33) 



Now i/^o is defined through = ^zJpO/ and i^o as a multiple of !/^0/ 



that fulfills 



!^0(±~)=±^. (35) 



We can then rewrite the expansion for the gauge field as 

A,, = ^+a^^;V+ + ?-3?'CV- + E4"Vn , ^z = o (36) 

n 

with 

1 

!/j±(z) = - ± i/^o such that i/'+(oo) =!/;_(— oo) = 1 and i/'+(—oo) = !/;_(oo) = (37) 



and 



fiH^'') = Pexp ^ - f^^ dz'A,{x^z') \ . (38) 



One stays in the Az = gauge under residual gauge transformations g{x^',z = 0) = h{x^'). Fixing the gauge to 
^_ = 1, we have the following gauge field expansion 16611 

A^{x,z) = U-\x)df,U{x)ip+{z)+l^BP{x)ip„{z), Az=0. (39) 

n>l 

The pion field is defined as 

U(xf') =Vexp\ - r dz'Az(xf',z')| (40) 



in order to use the same pion field as the one used in the sigma-model for low-energy effective QCD, (J(x'') = 
^2in{x>')/ fn ^ 7r(xf') = Tiat" , where the kinetic part of the action for the mesons is given by 

j d^x^ Tr(LI+a^!J)2. (41) 

Equating this action with the corresponding term in the Sakai-Sugimoto action after plugging in the expansion for 
the gauge field, 

^I'''iy ^ f = fj^'-^iy (42, 



leads to the identification 



f 2 ~ 3 ~ 3 / /■'^ V ^ ^ 



A ~ = 2 / dz— , (43) 

4 Mo "0 V -'0 {ul + uoz^Y'(> 



or 



A 3 / 'y^ 

ll(M..U„.) = -.M-i'-\ll ^^ („3^;,.).,J . (44, 
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where we have used the relations ijs}, ^ and the definition ( 11 , to determine the volume factor V in the action: 

V = ^kMI^^ (45) 

We now have all the ingredients to numerically fix the remaining parameters Mk, uq (thus L) and k. First, we 
determine k{Mk, mo) by demanding the constituent quark mass to be 0.310 GeV, 

m„{MK,uo,K) = Stz^mIk / du . = = 0.310 GeV =^ k(Mk, mq). (46) 

Ji/Mk n 1 

Then we use the experimental value for the pion decay constant to find uq{Mk), 

MMk,uo,k{Mk,uo)) = MMk,uo) = 0.093 GeV =^ mo(Mk). (47) 



Finally we solve the eigenvalue equation | [23) , which is now a function of Mk only, for m„^i. We refer the reader 
to the Appendix for more details. The value of Mk is then determined such that m„=i = nip = 0.776 GeV, the p 

meson being the lightest meson in the vector meson tower One identifies B^"~^^'' with p", {a = 1, 2), which can be 

recombined into the charged p^ mesons, Bj;""^''' with the neutral p^ meson, and bJ"-^^° with the cv^ meson. From 
these identifications it follows that mp = nico = nin^i in the Sakai-Sugimoto model. 

The results of our numerical analysis are 

Mk w 0.7209 GeV, ^ w 1.38 and k = w 0.006778, (48) 

Uk 2\bTC^ 

or 

Mk « 0.7209 GeV, L w 1.574 GeV"^ and k = w 0.006778, (49) 

21671^ 

where we used the formula |j3| describing the one-to-one relation between L and uq. The value found for L is 
approximately 2.8 times smaller than the maximum value of L, given by 

Lma. = ^ = « 4.358 GeV-1. (50) 
2 Mk 

From the values | |49| we do extract a relatively large 't Hooft coupling, A w 15. This allows us to ignore the Chern- 
Simons part of the action, which is a factor A smaller than the DBI-part. Nevertheless we will briefly comment 
on the contributions to the p meson mass equation originating from the Chem-Simons action in Section III A We 
also remark that for these values of the holographic parameters the numerical value for the effective string tension 
between a quark and an antiquark in this background, given by |68J 



Inoc 



V-^oo("Kkii(«K) «0.19GeV2, (51) 



is in good accordance with the value calculated on the lattice for pure SU(3) QCD, a k. 0.18-0.19 GeV^, as reported 
in Il84l 18511 . This is a nice illustration that the fixed values do have a reasonable QCD resemblance. 



C. Turning on a uniform magnetic field 

Under gauge transformations ^ e Lr(N^), the flavour gauge field transforms as 

A„(x'«)^gA™(x'")g-i+ga„,^'i (m = /(,z). (52) 



Since we have assumed the eigenfunctions ^nin > 1) to be normalizable (i/'n(z ±oo) = 0), the expansion |28l 
implicitly assumes we are working in the gauge A^, (z ±oo) = 0. The gauge potential can be made to vanish 
asymptotically by applying a gauge transformation g(x'',z) = U{x^,z), that cancels the asymptotic pure gauge 
configuration 

A,„{xf,z ±oo) = Ui^{x'-',z)dmU±{xf,z) (53) 
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that ensures a finite effective four-dimensional action. For arbitrary Njr > 2, the homotopy group for the functions 
(J : ]R4 U 00 ~ S4 ^> U(N|), — >• U{x^,z) is trivial, 7i/^{U{Nf)) = 0, so a continuously interpolating !J(N^)-valued 
function U{x^,z) that fulfills LJ(x'',z — > ±00) = U±{x^,z) can always be found. The case Nf = 2, which we 
consider, is an exception since 714(11(2)) = Z2. In the seminal paper of Sakai and Sugimoto Il66l , it was assumed 
that Nf 7^ 2. Yet it appears to be still possible to consider the gauge Am(z — > ±00) = 0. For Nf = 2, there will exist 

a 2 by 2 matrix function U{x^,z) interpolating between U+ and U- (if they are homotopic) or between U+ and U- 
(if U+ and U- are not homotopic), with (J_ defined as (73U-, so that U- is homotopic to the Z2 element =Fl if L?- 
is homotopic to the Z2 element ±1. The role of the multiplication is to switch sign of a row in U_, and thus also 
of the determinant. The sign of the determinant determines whether a U{2) matrix is homotopic to —1 or +1. Since 



the asymptotic pure gauge configuration 1 53 1 also equals U^^ {xf\ z)dniU±{xf' , z) , as can be easily verified, the gauge 
transformation g{x^',z) = U{xl^,z) will again cancel the asymptotic gauge potential. This argument'' extends the 
validity of the original Sakai-Sugimoto reasoning to the Njr = 2 case. 

One does not leave the gauge A,„ (z ±00) = under gauge transformations h G U{Nj:) that adopt x^-rndependent 
boundary values (/z+,/z_) = (limz^+co ^, lim^-s^-oo ^)- These boundary values of the residual gauge symmetry trans- 
formation h are interpreted as a global chiral symmetry transformation {hi,h^) E U{Njr)i x U{Njr)g^ in the dual 
QCD-like theory. By "lightly" gauging this chiral symmetry, i.e. making hi = = h{z ^ ±00) — h dependent on 
x'', one leaves the Am(z ±00) = gauge, the boundary value of the gauge field Am(z ±00) to be interpreted 
as an external background vector field A^, in the boundary field theory coupling to the quarks through a covariant 
derivative Vj, = dj, + A^, such that the Dirac action xpijfiVpiip remains invariant under local U{Nj:) transformations. 
Throughout this article we will be working in the Az = gauge instead, where we have the following expansion of 
the gauge field L67J (with our A^, equal to Ai = Ar = V in the notation of [67J) 

= (?+3,,?;i + ?+A,,?;i)!/;+ + (^-a^^I^ + ?_A,,Ci)i/;_ + E^^Vn- (54) 
Fixing the residual gauge symmetry within this gauge to (^^^ (x^) = ^_ (x''), the expansion becomes 

Af, =Af, + ^[TZ,df,Tz] + ^ {df,n - [n,Af,]) iPo + J^bI"^%. (55) 

In this paper we shall not look at the pions in detail, being mainly interested in the p mesons (identified with B^P), 
so essentially we will use the gauge field expansion 

=^/'+EsfVn- (56) 

To turn on an electromagnetic background field A™ in the boundary field theory we put (e being the electromagnetic 
coupling constant and Q^^ the electric charge matrix) 



A, = -ieQernA;:" = -ie j A™ = -ie Q12 + ^-a^j A^-, (57) 

which assigns the appropriate charge to the up- and down-quark. For the case of a constant external magnetic field 
along the X3-direction in the boundary field theory (F™ = diA™' = B), this amounts to setting 

A^ = -ieQemXiB6^i = -ix, J_^^ ^ 3^2 = + xteB3^2 (^-f) ' (58) 

A2 = xieB and A2 = A2/3 (59) 



or 



and 

where in the last line we defined the up- and down-components of the background field strength, Fj, and F^. 

^ We thank J. Van Doorsselaere for discussion on this point. 
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D. Effect of the magnetic field on the embedding of the probe branes 



We determine in this Section the eB-dependence of the embedding of the flavour D8-branes in the confining (we are 
working at zero temperature) D4-brane background Q. On each of the D8-branes lives an induced metric 



ds 



D8 



g^f,dx"'dx" {m,n = 0- 



/M\3/2 

\r) 



du' 



u^dCll 



or 



(„D8 „D8 „D8\ 
[800 'Sii 'Suu j 



fU\^n /'U\3/2 /MN3/2 



■■( 



RJ '\RJ 



1 / M N -3 

7 Vk. 



= (goo, Sii' Gun ) with Guu = gun + grr (3«t)^ 

and a gauge field A,,, for which we assume the background gauge field ansatz 



Afi — Afi = —ieQi,mXiB5^i2 (all other gauge field components zero). 



(61) 

(62) 
(63) 

(64) 



modeling an external magnetic field B = Be^ in the dual field theory. 
We plug the gauge field ansatz into the non-Abelian DBI-action 



Sum = -Ts d^x2 du €4 e'^ STr J- det [g^s + {2na')iFm„] 



(65) 



and solve for the embedding u' = du/ dr as a function of eB. To allow for the possibility that each of the two flavour 
D8-branes responds differently to the external magnetic field, we assume the following form of the metric in flavour 
space 







(66) 



the only difference between g^^ and g^^ being that the M-coordrnate appearing in follows the up-brane, varying 
from uq u to infinity, whereas the the w-coordinate appearing in g^^ follows the down-brane, varying from Mq^^ to 
infinity. This will turn out to generate a different embedding u' (eB) for up and down. We write 



,D8 



where we introduced the notation 



1' 



- "o,i() 

e{u - Mo,d) 



(67) 



(68) 



for our "generalized unity matrix" that indicates that everything multiplied by the Heaviside step function 
d{u — Uq u) (respectively 9[u — Uq ^)) will have to be integrated over u varying from Mo,m (respectively Uq ^j) to infinity 

The determinant in the action is 



det(^°^l' + i{2na')F) = det 
























gill' 


i{2na')Fn 















-i{2na')Fn 












X det S4 











^33!' 









V 













) 


Lorentz space 



detS4 X googng22g33Guu (l - (27Ta')^gn^g22^^i2) l' 

^ ^ 



detg X Au X e{u - Mo,u) 

detgxAdX0(M-Mo,d)y'fl„ 



(69) 



space 
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where we defined 



A = 



A„ 
A, 



(70) 



Since the matrix 1 69 1 is diagonal, the square root of it is equal to the matrix of the square roots of the diagonal 
elements and the STr reduces to an ordinary Tr, leading to the action 



Sdbi — -TsV^3+iV4gg ^ 



du u 



"0,u 



1 fu\-^ f nr- , 



f\R 



du u 



"04 



1 /M\-3 f 



f\R 



(71) 



with V]jj3+i = / d^x and u' as a function of eB to be determined for both the up- and down-brane. 
Omitting all the up- and down-indices for clearness, 

Sdbi ~ Sup + S^g^^„ 

/■CO 

S = -T8V]R3+iV4^7^ 2 / duu^ 

and using the short-hand 



1 /M\-3 / 



"0 



A 



+ fVA, 



V / 

we determine u' for each flavour from the conserved "Hamiltonian 



H = u' 



6u' 



, d^H = 0. 



(72) 



(73) 



(74) 



Expressing that this H is conserved and assuming a U-shaped embedding m' = at w = wq (with A(mo) and /(mq) 
denoted as Aq and /o): 



we find 



v//o 

,2 ^ 2 »V^ - »o/o^O 

" \r) ^ ulhA, ' 
reducing to the known U-shaped embedding for eB ^ whereby A ^ 1. 



(75) 



(76) 



1. Antipodal embedding (uq = Uf:^): no dependence on eB 
In the case uq — uk, we have /o = so 



0„.)a.(5)'i_^|AA. .0 (77, 

\uj f^u^fA-ulfoAo 

and the embedding function is constant, 

t(m) = t ~ 1, (78) 

(t = for the D8-branes and t = tt/ M^cl for the D8-branes, see the l.h.s. of Figure [l|, independent of the value 
of the magnetic field. In this case, there is thus no response of the chiral symmetry breaking to the magnetic field, 
a somewhat unphysical feature of the extremal Sakai-Sugunoto embedding, which is a direct consequence of the 
absence of a constituent quark mass in this setting. 
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2. Non-antipodal embedding (uq > Uf^): magnetic catalysis ofchiral symmetry breaking 
In the case mq > the embedding function for the flavour branes in the background is given by 

t(w) = T - 



Tu 
Trf 



with 



3uT 



R 



1 »o/o^o 

/2 w8/A - w8/oAo 



X 0(m — Uq), 



(79) 



(80) 



with Mg and A taking their up or down values in 3,(T„ or duT^i respectively. The up-brane and down-brane are thus 
no longer coincident in the presence of a magnetic field, as sketched in Figure |2j 




eB 



Uo,u(eB) 
Uo,d(eB) 




FIG. 2: The change in embedding of the flavour branes caused by the magnetic field eB models the chiral magnetic catalysis 
effect. The up-brane reacts the strongest to eB, corresponding to a stronger chiral magnetic catalysis for the up-quarks than for 
the down-quarks. 



The eB-dependent induced metric on the up- and down-brane is given by 



(81) 



and the action by 



S = -TsVj:,3+^Vig;^ 2 / du K3/2„5/2^^ 



with 



7b (w) 



u^fA-ulfoAo' 



(82) 



(83) 



We see that the non-Abelian DBI-action for the two D8-branes in the presence of a background magnetic field reduces 
to the sum of two Abelian actions. This represents the explicit breaking of the global chiral symmetry. 



U{2)l X U{2)r ^ {U{1)l X !J(1)r)" x (U(1)l x !J(1)r)'^, 
caused by the up- and down-quarks' different coupling to the magnetic field. 



(84) 
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The asymptotic separation L between D8- and D8-branes as a function of the magnetic field is 

L = 2 



:v//oAo/ rf^ , ^ ^ ^ (85) 



3 v"o ^0 



where we changed the integration variable to ^ = {u/uo)~^ ||47l, with j/j^ = uk/uq,!/ = u/uq and / = 1 — 



In Section II B the value of the geometric parameter L in zero magnetic field was determined at L(eB = 0) = 1.547 



GeV-i. We keep L = L(eB = 0) fixed while varying eB to determine Mq „/£((eB) and consequently, via 



, 9 /■"O 

ma{MK,uo,K) = 8Tz^Mj;^K du (86) 

J1/Mk 



(Mku)3 

the constituent quark masses nti, (eB) and mj{eB) of up- and down-quarks. 

In Figures [3] and |4] the numerically obtained dependence on eB of UQ^u/d' '^n ^rid are depicted. As uq rises with eB, 
the probe branes m the presence of the external magnetic field get more and more bent towards each other, driving 
them further and further away from the chirally invariant situation of straight branes, see Figure |2l This feature 
corresponds to a holographic modeling of the magnetic catalysis of chiral symmetry breaking [72] . As already 
mentioned, this "chiral magnetic catalysis effect" was already discussed for the Sakai-Sugimoto model in 1471 , albeit 
for a single flavour and without matching the free parameters onto QCD values. The constituent quark masses 



OTg(eB), which are related to the quantity Mo(eB) — m^/ accordingly increase (see Figure 
taking this chiral magnetic catalysis into account will translate into the p meson mass a' 
when ignoring the lowest Landau level shift (see next Section). 



if 

so 



, leading us to expect that 
so growing with eB, at least 



«o/«x 




FIG. 3: as a functiori of the magnetic field for the D8-brane corresponding to the up-quark, and the one corresponding to the 
down-quark. 



For small magnetic fields, we can compute the lowest order correction to mq(eB = 0) analytically, confirming the 
typical holographic (eB)^ dependence Il48l . To this end we approximate L by L = i'o(w) + (^S)^ ^i(m). The condition 
4(m) + (eB)2 li{u) = L{eB = 0) then has a solution of the form m = mq + (^B)^ mi with mq = 1-38 mk = 1-92 GeV"^ 

and Ml = — ^i(mo) (^^("o)^ ■ The corresponding constituent quark masses for small magnetic fields are 

mJeB) = 0.310 GeV + Sk^M\k / du- 



"0 Jl 



(AW 

1 
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We find 

OT„(eB) = 0.310 GeV + 0.582 (eB)2 + C(eB)3 and ;«d(eB) = 0.310 GeV + 0.145 (eB)^ + C(eB)3, (88) 

depicted in Figure |5] This quadratic dependence at small magnetic field is also encountered in other effective de- 
scriptions of the constituent quark mass, as in the PLSM^ model of [13 86 1, and although not explicitly mentioned, 
also an instanton based computation seems to give a quadratic-like power [19 1. Also the numerical data of [42 1 for 
the up- and down-quark chiral condensates are in accordance with a quadratic behaviour at small eB. We must how- 
ever mention that the latter lattice computations were done at nonvanishing current quark mass in an unquenched 
setting. Chiral perturbation theory predicts a linear behaviour [73] (see also the comments in [48J). Quenched lattice 
simulations of [37] confirmed this, although we notice that the small eB-behaviour does not seem to be precisely 
caught by the proposed linear fit. In fact, we are able to fit our result quite well with a linear fit if eB is not too large, 
see Figure |5] We used 

m„(eB)''"™'' w 0.303 GeV + 0.166 eB and mrf(eB)''"™'' w 0.301 GeV + 0.084 eB. (89) 

It is also instructive to see what happens at large magnetic field. As already pointed out in ||47| , we observe a 
saturation in Figure |4] Such a saturation was also seen for the first time using lattice simulations in ["42] for the up- 
and down-quark chiral condensates, carefully taken into account some subtleties related to an unphysical periodicity 
in the results, which is a typical lattice artefact. The results of ||T91 142| anyhow confirm the different response of 
the constituent up- and down-quark masses or chiral condensates to the magnetic field, a feature which we also 
reproduced here for the first time in the holographic Sakai-Sugimoto setting. It thus appears that our holographic 
results reproduce quite well the phenomenology of independent quenched QCD calculations. The unquenched top- 
of-the-bill simulations of |44| show a similar behaviour for the chiral condensate, at least at vanishing temperature. 
Similarly shaped curves for the chiral condensate extrapolated to the chiral limit can be found in |43|, which would 
be most relevant for comparison with our analysis, nevertheless keeping in mind that the lattice study |43J is for 2 
rather than 3 colours. 



III. EXPLORING THE p MESON MASS 

In this section we will determine the mass equation for the 5-dimensional gauge field fluctuation describing the 
four-dimensional p meson, again for both the antipodal and the non-antipodal embedding. The strategy is to plug 
the gauge field ansatz 

Ar = Ar + Ar (r = }l, u) 
^ T = T + f 

with (see l|64) and ||80)) 

A,, = -ieQemXiBSfij 



, , (91) 

Wo) 
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m, (GeV) nig (GeV) 




the same quantities. 



into the DBI-action describing the d5mamics of the flavour branes. We then expand the action to order {Ina')^ ^ 
and to second order in the fluctuations, and Anally integrate out the M-dependence to obtain an eflective four- 
dimensional theory for the p meson in the presence of the background magnetic field. Further we will use the 
notation t = T^tT"", f = fa", = A^^^t" and A,- = Aft", with all components field" real. 

We have to evaluate (for notational brevity we temporarily absorb the factor {2na') into the field strength) 

2 I du STr^-detK) =2 J du S^r^J - det{g^^ + iFrs), (92) 

with 

g?s^ = grs + grrDrTD.T, with D,- = 3, + [A,; ■] (93) 

and 

Frs=drAs-dsAr + [Ar,As]. (94) 

If we denote the argument a of the determinant (which runs over the Lorentz-rndices) as 

a = fl + fl(i) +flP) ^ 

with fl(") being n-th order in the fluctuations A, the determinant can be expanded to second order in the fluctuations 
as follows 

v/^detfl|^2 = ^/^dets|l + ^tr(rt(l)) + ^ (tr(fl-iflW))^- ^tr((fl-iflW)2) + ltr(fl-ifl(2))| ^ (95-, 
We denote the trace in Lorentz-space with a small tr, and the trace in flavour space with a capital (S)Tr. 

Splitting each component of a in its symmetric and antisymmetric parts 

a-^ = g + B 

ail) ^ ^(1) _^ (95) 

the expansion of the determinant < [95) to second order in the fluctuations can be written as 

V— detfl|^2 = V— detfl + V— detax 

{^tr(afl(i)) + 1 (tr(ga(i)))' - ltr{ga(^)ga(''> + Ba^^^'Ba^'^) + ^trl^a^^)) 

+ ^triBSiF) + ^ {tr{BStF)f - ^triGSiFgSiF + BS^FBSiF) + ^tr(S^2F) 

+ ^tr(afl(i))tr(^JiF) - ^tr(gfl(i'6JiF) - ^tr(g<5iFfifl(i' ) 1 . (97) 
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For our field ansatz we have 

ars = grs + gTT^i t3sT + iFr, 



4-^ = grr (a,T ([As,t] + Dsf) + ([A,,t] + D,f) d^r) , 



^iF,,=f(D,A,-D,A,)™'^"'"fF,, 



4? =gTT ([A„t] +D,f) ([As,t] +Dsf) +gTT ([A,,f]asT + 3,T[A„t]), 
The symmetric part Q of is diagonal. 



(98) 

(99) 

(100) 

(101) 
(102) 



Q = 



\ 



S33 



\ 



, with Guu = guu + gTT{duTf 



(103) 



(104) 



and the antisymmetric part B has non-zero components 

Bu = -B21 = i^iignSii^'^- 
The DBI-Lagrangian to second order in the fluctuations thus becomes 

STr e-'^^-deia\^2^^2^Af = Tr e^'^'V- deta + STr e-fg\^^JG^ugl^-/Ax 

r 1 1 

-G-^grr^uTfiig^^A-^ {{[Au,t] + D„f) F12 + ([Ai, t] + D^f) F2„ - ([A2, t] + D2f) Fi„) } 

where the notation for the factors g^^;}A~^|;/=i,2 coming from Q means that gj^^ is accompanied with a factor 
only for = 1,2. 

Approximating the action to order T^^ — (Ina')'^ ^ 1 /A^ while keeping in mind that every F carries a factor T^^ 
and that (see l|69)) 



lr2 



(105) 



A = l-g-2T-2F?2, 



(106) 



we have 



STr e-1'V-deta\^2^f2^Af,T-2 = Tr g-^V- detfl + STr e-'Pgl^^G^,glT~^^ 
{ " ^12) { ^ ([^.,t] + Duff + f] {GjgrrduT) 



-9^2 \ 1 



^=1 

^ -l„-lr2 l„-l^-lp2 



+ f 12^11 ^21 - f 12^11 {A\, A2] - F^,^ - ^g^^, G^u 

-GmlgrrduTFugn^ {{[Au,t] + D„t) F12 + ([Ai, t] + Dit) F2„ - ([A2, t] + D2f) Fi„) } . (107) 



A. Antipodal embedding (uq = uk): Landau levels : 



In the case uq = we have 



duf = 

T ^ 1 



Guu — 
[Ar,f] 



0, 



(108) 
(109) 



17 



simplifying | 107| drastically to a part in the scalar fluctuations and a part quadratic in the gauge fluctuations, given 

by 



-2x 



'DB7 



STr e-1'V-deta\^2j-2 = Tr e-'Pgl.^glj- 

[Piignhi - Fugn 1^1' ^2] - \g^lgvvflv - \g;,lgZuflu } , (110) 

where we also chose the gauge Au = 0. The DBI-action for the A", {a = 1, 2) components can then be written as^ 

2 

E 

/(,i/=i 

- \f2{plf{Mf - \h E F%esabP'"'p"V ] + Pion action (111) 
H,v=l J 

with/, (f = 1..3) the following functions of u: 

f,{u)=f,{u) = TsV^ilTza'fe-fgl^, (112) 
fiiu) = TsVii27za'fe-fglgngu^y\ (113) 




and where we retained only the first meson of the vector meson tower in the fluctuation expansion 1 55 1, A", = 
jO^(xf' )!/;(«) + pions {a = 1,2), because it is the most likely to condense first, being the lightest spin 1 particle, and 
J-'^y = Dj^pl — Dyp"^. There are no coupling terms between p mesons and pions at second order in the fluctuations 
in the DBI-action, which can be traced back to the different parity of ^o{z) = ipQ{u{z)) and ip{z) = ip{u{z)) (with 
m(z) = + U]^z^), ^o{z) being odd and i/'(z) even. In order to obtain a canonical kinetic term and mass term for the 
p mesons in the effective four-dimensional action, we impose that the ip{u) fulfill the standard conditions 

00 

dufixp^ = l, (114) 

> 

du f2{d,,^pf = ml. (115) 



These conditions combine to the eigenvalue equation |23i with J = f ^ for uq = Uf^. Per construction the lowest 



eigenvalue of this equation is = 0.776^ GeV^. The corresponding eigenf unction ip fulfilling the boundary condi- 
tions ip'{z = 0) =0 and ip{z — > ±00) = can be used to evaluate the last integral over u in the above action, which 
determines the non-minimal magnetic moment coupling k of the p mesons to the background magnetic field, 

du fyip^ = k, (116) 

UK 

related to the magnetic moment ji as ji = (1 + k)e/ (2m) so to the gyromagnetic ratio g as g — 1 -\-k ||87| . Because in 
this simple embedding we have /j = /s, we immediately see from the normalization condition that k = 1 and thus 
g = 2, describing a non-minimal coupling of the p mesons to the background magnetic field. 

The effective four-dimensional action thus takes the form of the standard four-dimensional action used to describe 
the coupling of charged vector mesons to an external magnetic field (i.e. the Proca action 1,87, ,88J or DSGS action for 
self -consistent p meson quantum electrodynamics to second order in the fields 1891 ): 



Seff 



I E l-li^^nf - Int^pip^ f - \ E fJvWV'") ■ (117) 



' Note that from here on contraction over Minkowski indices is assumed implicitly in the notation of squares, e.g. (f,",,)^ = fm,F,"i,7''''^ 
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This means that the Sakai-Sugimoto model with mq = automatically describes Landau levels for p mesons 
moving in an external magnetic field. Let us quickly repeat how to derive this from the effective action, following 
for example ||57| (up to conventions). 



The equations of motion for p^", 



combine to 



(118) 



(119) 



with T)pi = dj, + f = d^i + ieA''!^ for the charged combination p^, = (joj, + ip^ ) / -v/2, and the complex conjugate of 
this equation for the other charged combination ,0+ = (|D^, — ip^i) / Acting with D" on this equation of motion, 
and using the fact that [Dfj,Dv] = iV^^y, leads to the subsidiary condition 



DVi, = 0, 



which allows us to rewrite the equation as 



Vlipy - liFlypf - rn^py = 0. 



(120) 



(121) 



The condition | |120) and its conjugate are nothing else than the covariant (w.r.t. the electromagnetic background) 
generalizations of the usual Proca subsidiary conditions p^ = 0. 

Fourier transforming py e^^^'^~^^^ py, we find that the transverse combinations^ p^ ± ip2 respectively get a nega- 
tive or positive contribution ^eB to their effective mass squared as a consequence of their magnetic moment coupling 
to the magnetic field: 



E^{p^ ± ip^) = ik2 + xieBf +kl-dl + ml^ 2eB (p^ ± ip^ ) 



-d{ + {eBY xi 



eB 



+ kl + ml=f 2eB 



iPi ± iPi ) 



= eB{2N + l)+kl + ml^2eB {p- iip^). 
The following combinations of the transverse (w.r.t. B) polarizations of the charged spin-1 p meson fields. 



p = p^ + ip2 and p'^ = p^ 



(122) 



(123) 



have their spin aligned with the background magnetic field, S3 = 1, which decreases their energy with an amoimt 
eB. As a consequence, the energy or effective mass of the fields p and p^ in the lowest energy state, i.e. in the lowest 
Landau level N = and with zero momentum k^ along the direction of the magnetic field. 



eB, 



becomes imaginary when the magnetic field reaches the critical value of m^, 



eBc = OT^ w 0.60 GeV^, 



(124) 



(125) 



see Figure |6] We conclude from this that the field combinations ( 123| should experience a condensation, in 
accordance with ll30ll3T]| . 



' Because F^„, is only non-zero for ji,v 7^ 0, 3 there is no magnetic moment coupling term for the longitudinal components of the p fields, resulting 
in only the transverse components condensing according to the Landau equations of motion. 
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eB (GeV-) 



FIG. 6: The effective mass squared fJe^^ 
embedding of the flavour branes. ffJe^^ goes through zero at eBc- 



nip — eB of the field combinations p and p^ as a function of eB, for the case of antipodal 



Before concluding this Section, we mention that in principle there are also equations of motion for Af^, which, upon 
still ignoring the terms in the action, read (in the gauge A„ = 0) 



^0^ 



8„(Dy+)=0. 



(126) 



The latter equations are automatically fulfilled due to the Proca subsidiary conditions for the neutral mesons, the 
subsidiary condition 1 120) and its complex conjugate. 



B. Non-antipodal embedding (mq > mk): Chiral magnetic catalysis ^ eBc > nip 



Considering the case of non-antipodal embedding, uq > u^, corresponds to taking into account effects of the con- 
stituents of the mesons in the analysis, as the distance between uq and is related to the constituent quark mass 
through lEl. In contrast to the antipodal case, here the embedding of the flavour branes is affected by the magnetic 
field, as depicted in Figure [2] The eB-dependent change in embedding models two effects which both will have 
an influence on the mass ofThe p meson. Firstly, the rise of uq and consequently of the constituent quark masses 
ntu and with the magnetic field models magnetic catalysis of chiral symmetry breaking. As the constituents get 
heavier it is expected that the p meson will too, increasing the value of the critical magnetic field for condensation. 
Secondly, the fact that the up- and down brane are separated, models the explicit breaking of the global chiral sym- 

metry U{2) U{l)ii x /J(l)rf. The split between the branes generates another mass mechanism for the charged p 
mesons through a magnetically induced holographic Higgs mechanism: the charged ud and du combinations cor- 
respond to strings with one endpoint on the up-brane and the other endpoint on the down-brane, these will now 
stretch between the separated branes and because of their string tension will develop a (stringy) mass term. 
Taking into account the splitti ng o f the branes, t 7^ 1, severely complicates the analysis, mainly because the evalu- 
ation of the STr in the action 1 107 1 becomes quite technical in the case of non-coincident branes (the STr no longer 
reduces to a normal Tr). We are in the current paper moreover particularly interested in the effect of the chiral mag- 
netic catalysis on the p meson mass and the critical magnetic field. For this reason we will here approximate the 
embedding in Figure |2] by an embedding where the branes are still coincident in the presence of the magnetic field 
{duT 7^ but T ^ 1), joining at an averaged value uq of the holographic radius, determined from 



Laverage{uQ,eB) = I ( 85 ) with every A replaced by 



rage 



L(eB = 0) ^Mo(eB), (127) 



see Figure |7[ and |8] 

This form Or i-'average (mq/ eB) is the one you obtain when you postulate that the branes remain coincident in the pres- 
ence of the magnetic field by using = g^^l instead of | [67| |. The same approximation, i.e. not taking into account 
the splitting, is implicitly done when applying the magnetic field holographically by assigning a non-zero value only 

to A^, and not A^, (instead of (jsTj), as is quite often done in the literature, e.g. in [5¥-56|. In this approximation the 
Higgs effect on the p meson mass is ignored, but we will elaborate on this in a follow-up paper where we will handle 
the non-approximated case of splitted branes rigorously by suitably adapting tools developed in ll79ll8Tl . 
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FIG. 8: The average Mo(eB) (yellow) in the coincident branes approximation compared to uo ii{eB) (blue) and Mo_rf(eB) (red). 



In the case uq > and the current approximation of coincident branes we have 



X 9{u - Uq) 



P u»fA - ulfoAo 
f --^ 1 [Ay, t] = and Guu (^k t) drops out of the STr argument. 



(128) 
(129) 



simplifying 1 107 1 to a part in the scalar fluctuations and a part quadratic in the gauge fluctuations, given by 

STr e-1>V-deta\^2j-2 = Tr e^^v^^lj'^x 



(130) 



where we a gain chose the gauge A„ = and only retained the lowest mesons of the meson towers in the fluctuation 
expansions 1 55 1 for the gauge field and f (x^, u) =Y^„ U^"^ {x^)(pn («) for the scalar field: 



Al = plix^^'W) and f = aQ{x^),p{u) {a = 1,2). 



(131) 



Because both i/'(z) = ip{u{z)) and (p{z) = (p{u{z)) (with u{z) 



uqz ) are even functions, the last line of 1 107 1 



consists of integrals of the form J ^ dz{odd function of z} =0 and hence disappears. 

The DBI-action to second order in the charged gauge field fluctuations can then again be written as 



'DB7 



a=l 



2 

E 

//,l/=l 



d^x I du Y,\-\h{^;yf^'' - hiip^fiMf - \h E Fj.e3«bP''V"''/'4 +pion action (132) 
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with the functions /, (f — 1..3) dependent on both u and eB this time (through the eB-dependence of the embedding 
function 3„t in Gun)- 

fi{u,eB) =Mu,eB) = T^Viilna'fe-'Pgl^.jG^u, (133) 
h{u,eB) = T^V^{2na'fe-^g\gnGZ^i'^. (134) 

In order to obtain a canonical kinetic term and mass term for the p mesons in the effective four-dimensional action, 
we demand the !/'(«) to fulfill the standard conditions 



/ du fi{u,eB)ip^ = 1, (135) 

Ji/o(eB) 

rOO 

/ duf2{u,eB){du^pf = ml{eB). (136) 

JuQ(eB) 



These conditions combine to the eigenvalue equation 1 23 1 with 7 = "Y{eB) = ( 83 ) with every A replaced by Aaverage = 

2 

- m1/2^-i/23„(„5/2^-V29^^) ^ R^ml{eB)ip. (137) 

Numerically solving this eigenvalue equation for m^{eB) (see the Appendix for details) boils down to taking the ave- 
raged magnetic catalysis of chiral symmetry breaking into account, as we include the effect of eB on the embedding 
of the flavour branes, both through the induced metric on the branes and the changed value of Uq. The resulting p 
meson mass eigenvalue m^{eB) as a function of eB is depicted in Figure |9] As expected on grounds of Figure H it 
is an increasing function of magnetic field (where it should be noted that working in the approximation of comci- 
dent branes means that we consider the effect of an averaged increase for niqieB), equal for up and down, rather 
than the exact eB-dependences in Figure |4]|. The corresponding eigenfunction fulfilling the boundary conditions 
i/''(z = 0) = and ^{z ±00) = can be used to evaluate the last integral over u in the above action, which 
determines the gyromagnetic coupling constant k to be one, 

du f3{u,eB)ip^ = k = l, (138) 

UoieB) 

or g = A: + 1 = 2, again because of fi = /a- The effective four-dimensional action thus again takes the form of the 
standard four-dimensional action used to describe the coupling of charged vector mesons to an external magnetic 
field (i.e. the Proca action [87'. '885 or DSGS action for self -consistent p meson quantum electrodjrnamics to second 
order in the fields [.89J ), but with the influence of the constituent quarks reflected in the magnetic field dependence 
of rrir, 



Seff 



'^'^ E i-li^i'v)' - l^'pi^B)ip;f - \ E f JvWvj ■ (139) 



Completely analogous to the derivation in the preceding subsection, this action describes that the p and p^ fields in 
the lowest Landau energy level will get an effective mass 



E = m,ff= ^Jmj{eB) - eB, (140) 
which becomes imaginary when eBc = '«p(eBc), at 

eBc = m^ieBc) w 0.67 GeV^ w 1.1 m^, (141) 

see Figure [To] The averaged chiral magnetic catalysis pushes the critical magnetic field for condensation to a higher 
value, as expected. The increase is of the order of 10 percent. In future work, we shall probe the further increase in 
eBc due to the already alluded to Higgs effect caused by the flavour brane separation, to investigate to what extent 
we can get closer to the lattice estimate for the critical magnetic field, i.e. eBc of the order of 1 GeV^ ||40| . 
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FIG. 9: The p meson mass eigenvalue as a function of the magnetic field. 




FIG. 10: The effective mass squared "Jf^^ = mp(eB) — eB of the field combinations p and p^ as a function of eB, for the case of 
non-antipodal embedding of the flavour branes (blue), i.e. taking into accoimt the effect of chiral magnetic catalysis through the 
eB-dependenve of nip, compared to the case of antipodal embedding (red), m^jf'^ goes through zero at eBc- 



C. The Chern-Simons action and mixing with pions 

In principle, the DBI action ||5j needs to be complemented with a Chem-Simons piece which serves as the chiral 
anomaly in the QCD-like boundary theory: 

^-=2^^7^K^''-^^''+i^^')' ^'''^ 

where the notation implies wedge products of differential forms. Since it is a factor A smaller than the DBI-action 
and A w 15 ^ 1, we have ignored the CS-action in our above analysis. 

When one applies a magnetic field, one might intuitively guess that a quark and anti-quark forming a pion bound 
state might try to align their spins, thereby transforming into a p meson. In field theoretic terms, one might therefore 
expect a mixing between pions and p mesons due to a magnetic field. Given the different intrinsic parity of the 
participants in this mixing, this would necessarily correspond to an anomaly driven process, i.e. related to the Chern- 
Simons piece of the complete Sakai-Sugimoto action. It therefore looks interesting to investigate the CS-action to 
order meson gauge field squared in some more detail. We find that, after renormalizing by subtracting appropriate 
boundary terms 1.91 J , the action 

Scs ~ /Tr (e'""P'7'-A,„F„pF,, + 0{A^)) 

^ J 4Tr (A3F21F04 + A0F21F43 - A3F40F21 - A0F43F12 

+^2 [f 43^10 + ^10^43 + hlfos + ^03^41 + fsi^lO + ^40^31]) / 
does contribute interesting terms of the form pnB (see also |67|), 

~ B I |a[o/rV§] + I (a[o/r+P3] + d^o^^P^]) } , (143) 
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indeed describing couplings between p-mesons and pions (and other axial mesons that we have not taken into 
account). However, since there are no direct couplings between the transverse components {}i = 1,2) of the p fields 
and pions, we can conclude that, at the level of a OlA^) analysis, taking the CS-action into account will not affect 
the presence of the p meson condensate which is related to the transverse field components, see eq. | |123^ . 



IV. SUMMARY AND OUTLOOK 

We studied the p meson mass dependence on an external constant magnetic field, B = Be^, using the Nf = 2 
Sakai-Sugimoto model, which is a holographic dual of a QCD-like theory in the quenched approximation and chiral 
limit. We fixed all free parameters by matching them to known QCD values at zero magnetic field. As a first main 
result, we found that two non-antipodal embedded flavour branes do not respond equally strong to the magnetic 
field, which is the geometric translation of the up- and down-quark coupling to the magnetic field with different 
electromagnetic charges (respectively |e and — je), resulting in the up- and down-brane becoming separated at 
eB 7^ 0. This means that the chiral magnetic catalysis is stronger for the up-quarks than for the down-quarks. The 
presence of the magnetic field thus induces an explicit breaking of the chiral symmetry, which manifests itself in the 
non-Abelian DBI-action for the flavour branes reducing to the sum of two Abelian actions, one for the up- and one 
for the down-brane. As a result of the chiral magnetic catalysis, the constituent quarks gain mass (the up-quarks 
thus getting heavier than the down-quarks), and consequently also the p meson mass eigenvalue gets influenced. 

As a second main result, we have shown that the effective action for the p mesons in the presence of a magnetic 
background leads to Landau levels, whereby the effective masses of the combinations p^ =F ipf of transverse 
components of the p meson become tachyonic, indicative of a condensation of these charged fields. The value for 
the critical magnetic field at the onset of the condensation was estimated, taking into account the contribution 
to the effective p meson mass due to an averaged (i.e. equal for up and down) chiral magnetic catalysis, as 
eBc w 1.1 mj w 0.67 GeV^. Such a p meson condensation effect at very strong magnetic field was predicted and later 
confirmed on the lattice in the papers |30 31 40|. We will present the effect of the splitting of the branes on the p 
meson mass (i.e. the effect of the full non-averaged chiral magnetic catalysis plus a stringy Higgs effect) as well as 
the discussion of stability in the scalar meson sector in a follow-up paper ] 



Another step in this research program will be the explicit construction of the new vacuum in the condensed phase 
(either in the case of coincident or separated branes), to see to what extent the findings of ll30ll3Tll40ll can be further 
corroborated by holographic tools. Therefore, we shall need to discuss the condensation of p^ =F ip2 along the lines 
of [32 . 56J. It would also be interesting to include temperature effects, to study how the (super)conducting condensed 
p meson phase evolves. Hopefully, once these tasks are carried out, we shall be able to tell something more about 
a possible (T, eB) phase diagram of QCD as sketched in [92|. In principle, using the tools developed in this paper 
for the p mesons, also the axial vector mesons can be investigated. It would be interesting to find out whether these 
could condense or not, as they also couple to the magnetic field, both directly and indirectly. They were conjectured 
not to condense in li31J . Another line of investigation concerns the chiral restoration and deconfinement temperature 
structure of the Nf = 2 Sakai-Sugimoto model in the presence of the magnetic field. This will be presented elsewhere 

DMI. 

Appendix A: Solving the mass eigenvalue equation 

We numerically solve the eigenvalue equation | [23| with normalization condition | [2T) . Since 

j^-l/2y /2 _ y-1/2 fgj. „ _ 00, (Al) 

normalizability requires that 

j/7„(m = oo) = 0. (A2) 
For numerical purposes, it is easier to work on a compact interval. The following transformation 

m3 1 



Uq COS^ X 



(A3) 
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maps the interval under investigation onto xe [—n/2;n/2], with the boundary conditions ||A2} now reading 

i/^„(±7t/2) =0. (A4) 

The differential equation |[23) transforms into 



COS^ X 



t"0 — 
4 sm 



7\ 



„3 

-| COS^ X 



COS 



16/3 



■/o 



3x 



cos^''^ X 



sm 



„3 

-| COS^ X 



COS 



16/3 



■/o 



3x 



V 



i/'„(x) = R^\„Tp„{x), 



(A5) 



where we denoted the mass eigenvalue with A„. Due to the reflection symmetry x ^ — x we can split up the 
eigenf unction set in even/ odd ipn{x)'s and focus on the interval [0, n/l]. /Analogously as explained in 166, ,82 1, the 
even/ odd eigenfunctions correspond to odd/even parity mesons. We can thus demand that 



ip„{0)=0 or d^ip„{0)=0. 



(A6) 



For the p meson, we must look at the odd parity sector, in particular the even eigenfunction with the lowest eigen- 
value. We can temporarily replace ([23) and associated boundary conditions with the initial value problem 



Uq- 



COS^ X 



smx 



\ 



-| COS^ X 



COS 



16/3 . 



■/o 



dx 



COS^'''^ X 



sm 



-| COS"^ X 



COS 



16/3 



■/o 



dx 



V 



M3 



Aip,ip{0) =l,ip'{0) =0 (A7) 



where A is treated as a "shooting" parameter We numerically solved the previous differential equation for each 
value of A to give a unique ip\{x), consistent with the initial conditions. Since the coefficient functions appearing 
in | |A7) display a delicate behaviour --^ / around x = 0, some care is needed when using a numerical package, 
as this is precisely where we are imposing our initial conditions. The solution around x = can however be easily 
obtained using a Taylor expansion, and fed into the numerical procedure. Once the !/'a(^) is known, we can solve 
the equation 



'(A) = iPa{7t/2) 







(A8) 



for A, which is then precisely the mass eigenvalue of the original eigenvalue problem. 



For the eB-dependent eigenvalue problem ( |137| with normalization condition 1 135 1, the corresponding initial value 
problem can be written as I A7 1 with every /q — >■ /o^o / ^/ where in the case of the approximation of coincident branes 

A is replaced by Aaverage ~ — ^^-^ — i an d wh ere it is understood that each Uq appearing in the eigenvalue equa- 
tion is now eB-dependent, see UQ{eB) from (127 i. This can accordingly be solved numerically, using the explained 



shooting method, for the eB-dependent mass eigenvalue A(eB) = nip[eB). 
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